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Abstract Particle velocity distribution functions (VDF) in space plasmas often show non 
Maxwellian suprathermal tails decreasing as a power law of the velocity. Such distributions 
are well fitted by the so-called Kappa distribution. The presence of such distributions in 
different space plasmas suggests a universal mechanism for the creation of such suprather- 
mal tails. Different theories have been proposed and are recalled in this review paper. The 
suprathermal particles have important consequences concerning the acceleration and the 
temperature that are well evidenced by the kinetic approach where no closure requires the 
distributions to be nearly Maxwellians. Moreover, the presence of the suprathermal particles 
take an important role in the wave-particle interactions. 

Keywords Kappa distributions ■ space plasmas ■ kinetic models • waves and instabilities 



1 Introduction 

Nonthermal particle distributions are ubiquitous at high altitudes in the solar wind and 
many space plasmas, their presence being widely confirmed by spacecraft measurements 
dMontgomerv et al. ll968l ; lFeldman et al. ll97^|Pilipp et al. ll987l ; lMaksimovic et al. Il997al : 
Zouganelis 1 120081) . Such deviations from the Maxwellian distributions are also expected to 
exist in any low-density plasma in the Universe, where binary collisions of charges are suffi- 
ciently rare. The suprathermal populations are well described by the so-called Kappa (fc-) or 
ge neralized Lorentzi an velocity distributions functions (VDFs), as shown for the first time 
bv lVasvliunasl fl968). Such distributions have high energy tails deviated from a Maxwellian 
and decreasing as a power law in particle speed: 



V. Pierrard 

Belgian Institute for Space Aeronomy, av. circulaire 3, B-1180 Brussels, Belgium 

2 Center for Space Radiations, Universite Catholique de Louvain, chemin du cyclotron 2, B-1348 Louvain- 
La-Neuve, Belgium 

E-mail: viviane.pierrard@oma.be 

M. Lazar 

3 Research Department - Plasmas with Complex Interactions, 
Ruhr-Universitat Bochum, D-44780 Bochum 

E-mail: mlazar@tp4.rab.de 

4 Center for Plasma Astrophysics, Celestijnenlaan 200B, 3001 Leuven, Belgium 



2 



h [r ' v> iKiKwl^ r(K- i/2)r(3/2) >v ^ 2 Ki 

where vf K i = v (2fC — 3)&7}/ (Cm,- is the thermal velocity, m, the mass of the particles of 
species i, « their number density, J their equivalent temperature, v the velocity of the parti- 
cles, and r(x) is the Gamma function. The spectral index must take sufficiently large values 
K > 3/2 to keep away from the critical value K c = 3/2, where the distribution function {TJ 
collapses and the equivalent temperature is not defined. The value of the index K determines 
the slope of the energy spectrum of the suprathermal particles forming the tail of the VDF, as 
illustrated on Fig. 1. In the limit fc — > °°, the Kappa function degenerates into a Maxwellian. 
Note also that different mathematical definitions of Kappa distributions are commonly used 
and various authors characterize the power law nature of suprathermal tails in different ways. 

Here some comments are necessary because the conventional Kappa distribution func- 
tion m used to fit and describe high energy tails incorporates macroscopic parameters de- 
fined by the lo west moments o f the d istribution function. We follow some suggestive ex- 
planations from Hellberg et alj and first remember that w K t in (Q} was originally 
stated by I Vasvliunasl 1 1196a) to be the most probable particle speed. Thus, a characteristic 
(nonrelativistic) kinetic energy, W = mw 1 Ki /2, can be associated to the most probable speed, 
and by considering the second moment of the distribution function U = [d\ ffrnv 1 /2, the 
mean energy per particle reads W m = U /N = W(3k/2) /(fc — 3/2). Later. iFormisano et al. I 
dl973h have introduced a plasma temperature related to the mean energy per particle k^T = 
(2/3)W m = mwj.jK / (2k — 3), which is exactly th e equivalent temperature proposed bv lLeubneil 
d!982h and IChateau a nd Meye r- Vernetl d 1 99 lh by relating to the average energy kgT = 
m<v 2 >/3. On this basis, it was also shown that, in a Kappa distributed plasma, the Debye 
length is less than in a Maxwellian plasma, X K = X(2k — 3) /(2k — 1). Such a temperature 
definition, making use of equipartition of energy, although appropriate for an equilibrium 
Maxwellian distribution, is not strictly valid for a Kappa distribution, but there are practical 
advantages to using such an equivalent kinetic temperature, whic h can be a useful con cept 
already accepted in practice for non-Maxwellian distributions (see Hellberg et al. ( 2009) and 
references therein). 

Furthermore, generalizatio ns of thermodynamics based on th e Tsallis nonextensive en- 
tropy formalism dTsallisll995h have been used for several decades ( iLivadiotis and McComas I 
(2009) and references therein). The family of kappa distributions result from a new gen- 
eralized Lorentzian statistical mechanics formulated for a collisionless plasmas far from 
thermal (Boltzmann-Maxwell) equ ilibrium bu t containing fully developed turbulence in 
quasistationary equilibriu m ( Treu mannlll999allbl lLeubneJl2002t iFisk and Gloeckled 120061 : 
iGloeckler and F isk 2006; Treumann and Jaroschek 2008). Thus, in phase space, the power- 
law distributions of Kappa type describe marginally stable Gibbsian equilibria, and the pa- 
ramete r K controls the strength of the plasma particle correlatio n in the turbulent field fluctu- 
ations dHasegawa et al. |["l985l : lTreumann and Jaroschekll2008h . In such systems the temper- 
ature is redefined on the basis of a superadditive (superextensive) entropy because the inter- 
depen dence of subsystems contributes an extra amount to entropy (Treuma nn and Jaroschekl 
120081) . 

Such Kappa functions give the best fit to the observed velocity distribution functions, 
using only 3 parameters (the number density n, the temperature T and the parameter K char- 
acterizing the suprathermal tails). A sum of two Maxwellians can also represent distributions 
with enhanced suprathermal tails, but they need 4 parameters T\ and n-i, Ti representing 
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Fig. 1 The Kappa velocity distribution function for different values of the kappa parameter. 

Table 1 Comparison of different analytical expressions for a Maxwellian and a Kappa VDF. 
Parameter Maxwellian Kappa 

Number density n(r) = no exp ^— 
Temperature T(r) = Tq 

Escaping flux F(r) = "o"-(^'eJ w ) ex p ( 



«(r) = noW(l + f^)" + ' 

m = To^j- 2 (l + S) 
pi r \ _ "o-^t»'(l+"?/"' 2 ) 

V > 4(^-l)K 1 / 2 |l+r;/(Mr-)l' c 



the number density and temper atures of the two popul ations) and they generally give less 
good fits than Kappa functions (Zougane lis et al. 1120041) . 

Considering the suprathermal particles has important consequences for space plasmas. 
For instance, an isotropic Kappa distribution (instead of a Maxwellian) in a planetary or 
stellar exosphere leads to a number density n(r) decreasing as a power law (instead of expo- 
nentially) with the radial distance r and a temperature T increasing with the radial distance 
(instead of being constant), as shown by the expressions given in Table 1. R(r) is the poten- 
tial energy containing the effects of the gravitation, the electrostatic and centrifugal force. 
v e is the escape velocity and is the fraction of Gamma functions appearing in the Kappa 
VDF. Considering particles escaping as planetary or stellar wind, the Kappa distribution 
yield higher flux than a Maxwellian, since more suprathermal particles are able to escape. 

This review is organized in the following fashion. Reports of measurements or indirect 
detections of Kappa distributions in our interplanetary space are reviewed in the next section. 
In Section 3, we identify the mechanisms made responsible for the occurrence of nonthermal 
populations in different environments. Representative theories and scenarios developed for 
Kappa distributed plasmas are discussed in Section 4. In Section 5 we make a short overview 
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of the dynamics and dispersion properties of Kappa distributions including the recent results 
on the stability of anisotropic plasmas and kinetic instabilities. The impact and favorable 
perspectives for these distributions are discussed in Section 6. 



2 Detection of Kappa distributions 

Distributions with suprathermal tails have been observed in various space plasmas. Kappa 
distributions wit h 2 < K < 6 have been found to fit the observat ions and satellite data in 
the solar wind ijGloeckler et al. 19921: Maksi movie et al. 1 1l997al) . the terrestrial magne- 



tosphere dG loeckler and Hamilton 19S7). the terrestrial plasmasheet (Bameetal 
IChriston et al. 1 1 1 9881. 1 19891; lKletzing||2003h. the magnetosheath (jFormisano et al 
dPi 



the radiation belts ( Pierrard and Lemaire 



1967; 
19731) . 



1996bl; IXiao eUuH IIoOSc). t he magnetosphere of 



other planets like Mercury dChristonll987 ), the plasmasheet of Jupiter ([Collier and Hamilton I 

19951), the magnetosphere of Jupiter dKrimigis et al. ll98ll : lMauk et al. 12004)), Sat urn (Krimigis et al. I 
198a-|Schippers et al. l2008llDialvnas et al II2009I). Uranus dKrimigis et al. fl986h . Neptune 



1 Mauk et al. 1 199 lb , and even on Titan jDe la Have et al. l2007t) and in the Io plasma torus as 



observed by Ulysses (Mever-Vernet et al. 1995), Cassini ( Steffi et al. 2004) and the Hubble 
Space Telescope ( Retherfor d et al. 112003b . 

In the solar wind, electron velo city distributions are c haracterized by a thermal core 
and a halo suprathermal population dPierrard et al. 112001a) , as illustrated on Fig. [2] These 
electron VDF are also characterized by a strahl component aligned with the interplanetary 
mag netic field. Electron VDF measured by Ulysses have been fitted with Kappa functions 
bv iMaksimovic et al. I dl997ah . They show a global anticorrelation between the solar wind 
bulk speed and the value of the parameter K, that supports the k inetic theoretical result tha t 
the suprathermal electrons influence the solar wind acceleration dMaksimovic et al.l[l99 7b). 
Solar wind parti cle VDF observed b y CLUSTER have also been fitted by the generalized 
Kappa function dOureshi et al. II2002I) . Radial evolution of nonthermal electron populations 
in the low-latitude solar wind with Helios, Cluster, and Ulysses observations shows that 
the relative number of strahl electrons is dec reasing with radial distance, whereas the rela- 
tive number of halo electrons is increasing dStverak et al. 1 12009). Observations of electron 
suprathermal tails in the solar wind suggest their existence in the solar corona, since the 
electron mean free path in the solar wind is around 1 AU. The ion charge measurements 
stated by Ulysses were found to be co nsistent with co ronal Kappa VDF of electrons with 
kappa index ranging between 5 and 10 dKo et al. II 19961) . 

To be able to measure the suprathermal electron parameters in space plasmas, the quasi- 
thermal noise spectroscopy was implemented with Kappa distrib utions using in situ Ulysses/ 
URAP radio measurements in the solar wind dZouganelis Il2008h . This noise is produced by 
the quasi-thermal fluctuations of the electrons and by the Doppler-shifted thermal fluctua- 
tions of the ions. A sum of two Maxwellians has extensively been used but t he observations 
have shown that the electrons are better fitted by a kappa distribution function dLe Chat et al. I 
l2Q09h . 

Solar wind ion ( 20 Ne, 16 and 4 He) distribution functions measured by WIND and aver- 
aged over several days have also been fitted by Kappa functions dColrieretal.il 19961) . Low 
values of K (between 2.4 and 4.7) are obtained due to the presence of high suprathermal 
tails. Suprathermal solar wind particles were also measured in H + , He ++ , and He + distri- 
buti on functions d uring corotating interaction region (CIR) events observed by WIND at 1 
AU dChotoo et al. Il2000l) . 
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Fig. 2 Velocity distribution function observed by the satellite WIND at 1 UA (215 Rs) for high speed solar 
wind electrons. Bottom panel: Isocontours in the plane of velocities (normalized to the thermal velocity) par- 
allel and perpendicluar to the interplanetary magnetic field. Top panel: Parallel (solid line) and perpendicular 
(dotted line) cross section of the observed VDF. The dashed line represents the Maxwellian distribution that 
well fits the core of observed VDF. jPierrard et al. 119991) 



3 Generation of Kappa distributions in space plasmas 

Various mechanisms have been proposed to explain the origin of the suprathermal tails of 
the VDFs and occurrence of Kappa-like distrib utions in the solar wind, the corona and other 
space plasmas. The first one was suggested by Hasega wa et al. I dl985h who showed that a 
plasma immersed in a suprathermal radiation field suffers velocity-space diffusion which is 
enhanced by the photon-induced Coulomb-field fluctuations. This enhanced diffusion uni- 
versally produces a power-law distribution. 

ICollierNl993h uses random walk jumps in velocity space whose path lengths are gov- 
erned by a power or Levy flight probability distribution to generate Kappa-like distribution 
functions. T he adiaba t ic tran sport of suprathermal distributions modeled by Kappa functions 
is studied in Collier ( 19951) . The same author shows that space plasmas are dynamic sys- 
tems where th e energy is not fixed, so that the maximum entropy should not be considered 
dCollier 112004b 

iTreumannl d200lh developed a kinetic theory to show that Kappa-like VDFs correspon d 
to a particular thermodynamic equilibrium state dTreumannl 1 1 999bl : iTreumarm e.t al .1120041) . 
A new kinetic theory Boltzmann-like collision term including correlations was proposed. 
In equilibrium (turbulent but stable state far from thermal equilibrium), it yields the one- 
particle distribution function in the form of a generalized L orentzian resembling but not 
being identical with the Kappa distribution jTreumannlll999ah . 

iLeubnerl J2002h shows that Kappa-like distributions can result as a consequence of the 
entropy generalization in nonextensive Tsallis statistics dTsallislll995h . physic ally related to 
the long range nature of the Coulomb potential, turbulence and intermittency dLeubner and Vorosl 
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l2005l :l Treumann and Jaroschekll2008h . The Kappa distribution is equivalent to the q distri- 
bution function obtained from the maximisation of the Tsallis entropy. Systems subject to 
long-range interactions a nd correlations are fundamentally related to non-Maxwellian dis- 
tributions dLeubner| [2008). Core -halo distributio n functions are a natural equilibrium state 
in generalized thermostatistics ( iLeubnej [2004b ). Fundamental issues on Kappa distribu- 
tions in space plasmas and interplanetary proton distributions are emphasized in iLeubnerl 
l l2004al) . rLivadiotis and McComas I d2009l) also examined how Kappa distributions arise nat- 
urally from Tsallis statistical mechanics and provide a solid theoretical basis for describing 
complex systems. 

The generation of suprathermal electrons by resonant interaction with whistler waves i n 
the solar corona and wind was sugg ested bv lVocks and Mann I d2003l) ; IVocks et al. I ( l200i) . 
Introducing antisunward-propagating whistler waves into a kinetic model in order to provide 
diffusion, their results show that the whistler waves are capable of influencing the solar wind 
electron VDFs significantly, leading to the formatio n of both the halo a nd strahl populations 
and a more isotropic distribution at higher energies dVocksetal.ll2005l) . 

In an ambient quasi-static magnetic field, plasma charges gain energy through the cy- 
clotron resonance and the transit time damping (magnetic Landau resonance) of the linear 
waves. This is the cas e of high frequency whis tler mode that enhance the energy of electrons 
in Earth's foreshock dMa and Summ ers 1998), or that of MHD waves wh ich can accelerate 
both the electrons and the protons in the solar flares ( Mille d 1991L 1 19971) . and in the inner 
magnetosphere dSummers and Ma 2000). 

When large amplitude waves are present, t he nonlinear Landau damping can be respon - 
sible for the energization of plasma particles dMillerill99ll : lLeubneJ[200bIlShizgalll2007h . 
Stochastic acceleration of plasma particles in compressional turbulence seems to be consis- 
tent with the power law spectra occ ured throughout the helio s heath downstream from the 
termination shock of the solar wind dFisk and Gloeckler 20061: iFisk and Glocklejl2007l) . A 
mechanism for the generation of electron distribution function with suprathermal tails in the 
uppe r regions of the solar atmosphere in the presence of collisional damping was suggested 
bv lVinas et al. I d2000h as due to finite-amplitude, low-frequency, obliquely propagating elec- 
tromagnetic waves. The nonthermal features of the VDFs can also result from superdiffu- 
i it — ; — -i 

sion processes (Treumann 1997), and due to heat flows or the presence of the temperature 
anisotropies dLeubner and Vinaslll986T) . 

In the same spirit, iMa and Summers! < fl99^> considered the steady state solution of the 
Fokker-Planck (FP) equation and obtained a Kappa distribution for a quasi-linear wave- 
particle diffusion coeffici ent that v aries i nversely with the particle speed for velocities larger 
than the thermal speed. IShizgall d2007t> used the same FP equation to study the relative 
strengths of the wave-particle interactions and Coulomb collisions. The formation of high- 
energy tails in the electron VDF was also investigated with a FP model bv iLie-Svendsen et al. I 
dl997l) . 

Note that a one dimensional, electrostatic Vlasov model has been proposed for the gener - 
ation of suprathermal electron tails in solar wind conditions dCalifano and Maneenev 120 08). 
The po ssible development of Kappa velocity distribution was also illustrated in lHau and Fu I 
d2007l) by the problem of low-frequency waves and instabilities in uniform magnetized plas- 
mas with bi-Maxwellian distribution. 

Whatever the mechanisms of suprathermal tails formation, the kappa function is a useful 
mathematical tool to generalize the velocity distributions to the observed power law func- 
tions, the particular Maxwellian VDF corresponding to the specific value of K — > °°. 
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4 Theories based on the existence of Kappa distributions 



4.1 Star's corona 



IScudderl ( tl992iJ ibh was pioneer to show the consequences of a postulated nonthermal distri- 
bution in stellar atmospheres and especially the effect of the velocity filtration: the ratio of 
suprathermal particles over thermal ones increases as a function of altitude in an attraction 
field. The anticorrelation between the temperature and the density of the plasma leads to 
this natural explanation of velocity filtration for the heating o f the corona, without deposit- 
ing wave or magnetic field energy. Scudder dScuddei 1992b) determined also the value of 



the kappa parameter for different groups of stars. IScuddei (1994) showed that the excess of 



Doppler line widths can also be a consequence of non thermal distributions of absorbers and 
emitters. The excess brightness of the hotter li nes can satisfactorily be accounted for by a 
two-Maxwellian electron distribution function dRalchenko et al. 1120071) and should be also 
by a Kappa. Note that many solar observations implicitely assume that the velocity distribu- 
tions are Maxwellian in their proper frame, so that the presence of suprathermal tails should 
lead to reinterpretation of many observations. 

The mechanism of velocity filtration in solar corona ha s been propose d to explain the 
high energy electrons at higher altitudes in the solar wind dScuddedll992l lbl). Velocity fil- 
tration was also applied to heavy ions in the corona to explain the ir temperatures more tha n 
proportional to their mass observed in the high speed solar wind dPierrard and Lamv l2003h . 

Studying the heat flow carried by Kappa distributions in the solar corona. lporelli and Scudderl 
d 19991) demonstrated that a weak power law tail in the electron VDF can al low heat to flow up 
a radia lly directed temperature gradient. This result was also confirmed bv lLandi and Pantellini I 
d200lh who obtained the heat flux versus fc in a slab of the solar corona from a kinetic simu- 
lation taking collisions into account. For K > 5, the flux is close to the Spitzer-Harm classical 
collisional values while for smaller values of fc, the heat flux strongly incr eases and changes 
of sig n! If K is small enough, the fast wind can be suprafhermally driven dZouganelis et al. I 
2005). This shows the inadequacy of the classical heat conduction law in space plasmas 
and the importance to deal with non-Maxwellian velocity distribution such as Kappa VDF 
( Mever- Vernet 1ll999l , l2007l) . 

Note that the Kappa distribution is also consist ent with mean electron spectra produc- 
ing hard X-ray emission in some coronal sources l lKasparova and Karlickv ll2009h . More- 
over, the low coronal electron temperatures and high ion charge states can be reconciled if 
the coronal electron distribution function sta rts to develop a significant suprathermal halo 
already below 3RS dEsser and Edgar 1 12000). Effects of a Kappa distributi on function of 
electrons on incoherent scatter spectra were studied by ISaito et aL ] d2000h . The equilib- 
rium ionization fractions of N, O, Ne, Mg, S, Si, Ar, Ca, Fe and Ni were calculated for 
Maxwellian and Kapp a VDF based on a balance of ionization and recombination processes 
dWannawichian et al. 1120031) for typical temperatures in astrophysical plasmas. Low kappa 
values lead generally to a higher mean charge. 

The Coulomb focusing effects on the bre msstrahlung spec trum are investigated in aniso- 
tropic bi-Lorentzian distribution plasmas in lKim et al. I d2004l) . Plasma screening effects on 
elastic el ectron-ion collis i on pro cesses in a Lorentzian (Kappa)-distribution plasma are an- 
alyzed in ljung and Hong l d2000h . 
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4.2 Solar wind 

IPierrard and Lemaiie ( 1996a) developed a kinetic model of the ion - exosp heres based on the 
Kappa VDF. The heat flux was sp ecified inlPierrard and Lemaire I dl998h . Their model has 
been applied to the solar wind by iMaksimovic et al. I I I 1997a) and predicts the high speed 
solar wind velocities with reasonable temperatures in the corona and without additional ac- 
celeration mechanism. The presence of suprathermal electrons increases the electrostatic 
potential difference between the solar corona and the interplanetary space and accelerates 
the solar wind. The collisionless or weakly collis ional models in corona jScuddej[l992bl : 
IMaksimovic et al. lll997bl : IZouganelis et al. l2005h . all using VDFs with a suprathermal tail, 
are able to reproduce the high speed streams of the fast solar wind emitted out of coro- 
nal regions where the plasma temperature is smaller, as well as the low speed solar wind 
originating in the hotter equatorial regions of the solar corona. 

The exospheric Lor entzian (or Kappa) m odel was extended to non monotonic potential 
energy for the protons jLamvetal.ll2003cJ) and shows that the acceleration is especially 
large when it takes place at low radial distances in the coronal holes where the number 
densit y is lower than in other regions of the corona, as illustrated on Fig. [3] Zouganelis et al. 
(2004) demonstrated with a parametric study that this acceleration is a robust result produced 
by the presence of a sufficient number of suprathermal electrons and is valid also for other 
VDF with suprathermal tails than Kappa. 

The acceleration of the solar wind heavy ions is investigated in IPierrard et all ( 120041) : 
due to their different masses and charges, the minor ions reach different velocities. Even if 
their mass on charge ratio is always larger than that of the protons, they can be accelerated 
to velocities larger than that of the protons if their temperatures are sufficiently high in the 
corona. 

Adding the effects of the Coulomb collisions, a kin etic solar wind model based on the 
solution of the Fokker-Planck equation was developed ( IPierrard et al. lfl999l l2001ah . Typi- 
cal electron VDFs measured at 1 AU by WIND have been used as boundary condition to 
determine the VDFs at lower altitudes and it was proved that, fo r several solar radius , the 
suprathermal populations must be present in the corona as well ( IPierrard et al. II 19991) . In- 
deed, since the particle free path increases as v in a plasma due to the properties of Coulomb 
collisions, the suprathermal particles are non collisional even when thermal particles are sub- 
mitted to collisions. High energy tails can dev elop for Knu d sen numbers (i.e. ratio of mean 
free path over scale height) as low as 10 3 (IShoublll983l) . lMarsch and Livi 1 dl985l) have 
studied the collisional relaxation process and the associated rates (diffusio n and friction) 
for a n onthermal solar wind with Kappa VDFs. The Fokker-Planck model ( IPierrard et al. I 
2001a) illustrated also the transformation of the VDF of the electrons in the transition re- 
gion between the collision-dominated region in the corona and the collisionless region at 
larger radial distances. The VDF became more and more anisotropic in the transition re- 
gion. Contrary to exospheric models that are analytic, these collisional kinetic models have 
to be solved numerically. The spectral method of expansion of the solution in orthogonal 
polynomials converges faster for sup rathermal plasmas by using p olynomials based on the 
Kappa function weight developed by Magnus a nd Pierrard I d2008h . 

As a result of low collision rates in the interplanetary plasma, the electrons and the ions 
develop temperature anisotropies and t heir VDFs become skewed and develop tails and heat 
fluxe s along the ambient m agnetic field (Mars ch et al. 1 19 82: PiliDD et a l. ll987l : ISalem et al. I 
l2003l :l Stverak et al. 120081) . Moreover, field-aligned fl uxes of (s uprathermal) particles can be 
encountered at any altitude in the solar wind JPihppet all 19901) . but they become prominent 
in energetic shocks, like the coronal mass ejections or the fast solar wind at the planetary bow 
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Fig. 3 Number density, electrostatic potential, bulk speed and potential of protons for different solutions 
of the Kappa kinetic model of the solar wind. The lowest blue velocity curve (u=320 km/s at 215 Rs) in the 
bottom left panel corresponds to a Maxwellian model with a starting radial distance (called exobase) at tq = 6 
Rs. The middel red line (u=460 km/s at 215 Rs) corresponds to a Kappa VDF with K = 3 and tq = 6 Rs. The 
upper black line ( u=640 km/s a t 215 R s) corresponds to a Kappa VDF for K = 3 with an exobase at ro = 1.2 
Rs (Adapted from lLamv et al. I <2003bl) ). 



shock, giving rise to counterstreaming plasma events dFeldman et al. Hl974l : lGoslingetal.l 
1 19931 : ISteinberg et al~1l2005l : iGloeckler and Fiskll2006h . The same mechanisms mentioned 
above can induce anisotropic nonthermal VDFs in counterstreaming plasmas, and such com- 
plex plasmas can hold an important amount of free energy that makes them unstable to the 
excitation of waves and instabilities. 

iFainberg et al. I dl996h suggested that the non-thermal electrons can contribute as much 
as 50 % of the total electron pressure within magnetic clouds. Using this hvpothesis. lNieves-Chinchilla and Vinas I 
d2008h have assumed two populations of electrons in magnetic clouds: a core Maxwellian 
and a halo Kappa-like distribution. They found that kappa values exhibit either minor dif- 
ferences or, for some events, can be greater inside than outside the magnetic clouds. 

Observations from VOYAGE R indicate that ions in the outer heliosphere are well de- 
scribed by Kappa functions dPecker et al. 1 120051. The effects of a Kappa distribution in the 
helioshea th on the global hel i osphe re and energetic neutral atoms (ENA) flux have been 
studied in Heerikhui sen et al. I d200St) . The use of Kappa, as opposed to a Maxwellian, gives 
rise to greatly increased ENA fluxes above 1 keV, while medium energy fluxes are somewhat 
reduced. The effect of a Kappa distribution on the global interaction between the solar wind 
and the local interstellar medium (LISM) is generally an increase in energy transport from 
the heliosphere into the LISM, due to the modified profile of ENA's energies. This results in 
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a motion of the termination shock (by 4 AU), of the heliopause (by 9 AU) and of the bow 
shock (25 AU) farther out, in the nose direction. 



4.3 Earth's exosphere 

The Kappa model of ion-exosphere dPierrard and Lemaire Ill996bh has been used to study 
differ ent plasma regions in the magnetosphere of the Earth. Pierrard dPierrard 1 1996l : |Pierrard et al~l 
l2007h obtained new current- voltage relationships in auroral regions when suprathermal par- 
ticles are assumed to be present with Lorentzian and Bi-Lorentzian distributions. Field- 
aligned conductance values were also estimated from Maxwellian and Kappa distr ibutions 
in quiet and disturbed events using Freja electron data dOlsson and Janhunen II 1998h . 

Introducing a Kappa model appears to resolve discrepancy between calculations and 
observations of resonant plasma echoes and emissions used for in-situ measuring the lo- 
cal electron densi ty and the magnetic field strength in the magnetospheric environments 
dVifias et alj|2005l) . 

The three dimensional plas masphere has been modele d using Kappa velocity distribu- 
tion functions for the particles llPierrard and Stegen H 2008): this physical dynamic model of 
the plasmasphere gives the position of the plasmapause and the number density of the parti- 
cles inside and outside the plasmasphere. The effects of suprathermal particles on the tem- 
peratu re in the terrestrial plasmasphere were illustrated using Kappa functions in lPierrard and Lemaire I 
il200lh . 

The terrestrial polar wind is in some way similar to the escape of the solar wind: sim- 
ilar effects of suprathermal particles appear a nd lead to an increase of the escaping flux 
dLemaire and Pierrard 11200 H ; iTam et al. Il2007h . Along open magnetic field lines, the wind 
speed is inc reased by the presence o f suprathermal particles. A Monte Carlo simulation de- 
veloped by Bargho uthi et al. I d200lh shows the transformation of H + polar wind velocity 
distributions with Kappa suprathermal tails in the collisional transition region. 



4.4 Planetary exospheres 

Mever-Vernet 1 d200lh emphasized the importance of not being Maxwellian for the large 
structure of planetary environments. For bound structures shaped along magnetic field lines, 
the temperature increases with the distance, in contrast to classical isothermal equilibrium 
(see Table 1). The rise in temperature as the density falls is a generi c pr operty of distribu- 



tions w ith suprathermal tails, as shown in lMever-Vernet et al. 1 1 1995h and lMoncuquet et al. I 
(2002) to explain the temperature inversion in the Io torus. 



The ion-exosphere Kappa model ( Pierrard and Lemaire 1996a) has been adapted to the 



Saturnian plasmasphere (Moore and Mendillo I I2OO5I) . The kappa index gives an additional 



parameter to fit observations of Cassini. The polar wind and plasmasph ere of Jupiter and 
Saturn were also recently modeled with Kappa functions dPierrard 1 12009): the suprathermal 
particles increase significantly the escape flux from these giant planets, so that the iono- 
sphere become an important source for their inner magnetosphere. 

Spacecraf t charging environments at the Earth, Jupiter and Saturn were also obtained by 
Garrett and Hoffman ( 2000) using Kappa distributions for the warm electrons and protons. 
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5 Dispersion properties and stability of Kappa distributions 

In many circumstances, the wave-particle interactions can be made responsible for estab- 
lishing non-Maxwellian particle distribution functions with an enhanced high energy tail 
and shoulder in the profile of the distribution function. In turn, the general plasma dynamics 
and dispersion properties are also altered by the presence of nonthermal populations. Thus, 
the waves and instabilities in Kappa distributed plasmas, where collisions are sufficiently 
rare, are investigated using kinetic approaches based on Vlasov-Maxwell equations. 



5.1 Vlasov-Maxwell kinetics. Dielectric tensor 

Using a kinetic approach. ISummers etal.ldl994l) have calculated the dielectric tensor for the 
linear waves propagating at an arbitrary angle to a uniform magnetic field in a hot plasma 
with particles modeled by a Kappa distribution function. Despite the fact that the elements 
of this tensor take complicate integral forms, this paper is of reference for the theory of 
waves in Kappa-distributed plasmas. This dielectric tensor can be applied for analyzing the 
plasma modes as well as the kinetic instabilities in a very general context, limited only by the 
assumptions of linear plasm a theory. The analytical dispersion relations derived previously 
( iThorne and Summer sl l 19861) in terms of the modified Bessel functions of the lowest order 
/o, I\, Kq, and (which are tabulated) were restricted to a weak damping or growth of 
plasma waves by resonant inte ractions with plasma p articles. 

The approach developed in lSummers et al.l d 1994) also restricts to the distributions func- 
tions which are even functions of parallel velocity of particles, vii (where parallel or perpen- 
dicular directions are taken with respect to the stationary magnetic field), and this is what 
the author s called an usual cond ition in practice. This condition fails only in some extreme 
situations dSummers et al.lll994h as are, for instance, the asymmetric beam-plasma struc- 
tures developing in astrophysical jets or more violent shocks. These cases can however be 
approached distinctivel y in the limits of waves propagating along or perp endicular to the 
ambient magnetic field dLazar et al. ll2008al . l2009al : lLazar and Poedts lEoTM) . 



5.2 The modified plasma dispersion function 

Early dispersion studies have indeed described the simple unmagnetized plasma modes 
and the field aligned waves in the presence of an ambient magnetic field dLeubnej[l983l ; 
ISummers and Thornelll99ll) . Let us first remember the most important analytical changes 
introduced by the power law distributions of plasma particles. 

Kinetic theory of an equilibrium Maxwellian plasma naturally produce a dispersion ap- 
proac h based upon the well-know Fried and Conte plasma dispersion function dFried and Contel 

Z(/) = ir^^, Im(/)>0. (2) 



n l l 2 J-<» x-f ' J kyikBT' 

For a n on-Maxwellian plasma characterized by the Kappa distribution function dTb. lSummers and Thornd 
dl99lh have derived a modified plasma dispersion function 
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where the spectral index fc was first restricted to positive integers fc > 1.5 . This new disper- 
sion fu nction has been ingenuously generalized to arbitrary real fc > 1 .5 by Mace and Hellberg 
dl995h also showing that it is proportional to the Gauss hypergeometric fun ction, 2F1 . Exten- 
sive characterizations for th i s Kap p a dispersion function have been made by Summers and Thome 

199l h : lMace and Hellberel dl995h : lMacld2003h : IValentini and D'Agostall2()07h : lMace and H ellberg 



20091) . Furthermore, for a (space) plasma immersed in a stationary magnetic field, which 



determines a preferred d irection for the acceleration and motion of electrons and ions, 
iHellberg and Macel d2002l) have introduced an hybrid Kappa-Maxwellian distribution func- 
tion, one-dimensional Kappa along the magnetic field lines and Maxwellian perpendicular 
to this direction. The new dispersion function Z k m obtained for this anisotropic plasma dif- 
fers fromZK in <[3j> in the power to which the term (1 +x 2 /k)~ l is raised: in the isotropic case 
it is K + 1 , while in the o ne-dimensional case it is K. Relation between these two integral 
functions takes two forms dHellberg and Macell2002l :l Laza r et al. I l2008bh 



(y-1) 3 / 2 
(^-3/2)^/2 

f 2 



Zr-1 



K-1 



1/2 



i + -)z,(/) + ^ 



1 

2k 



(4) 



and both Z k m{/) and Z K (f) functions approach the Maxwellian dispersion function Z(/) 
from lf2} in the limit of fc — > °°. 



5.3 Isotropic Kappa distributions 

The effect of an isotropic Kappa popul ation on plasma modes has been described by the first 
dispersion studies already reviewed bv IHellberg et~aT] d2000bh and lHellberg et al. I d2005h . 

5.3.1 Unmagnetized plasma 

The fluctuations are, in general, enhanced in low-fc plasmas, and the electromagne tic and 

electrostatic dispersion relations show significant dependence on the spectral in dex K dThorne and SummersI 
ll99ll ; ISummers and Thornell 1 9921 ; iMace and Hellberglll995r . lMace et al.lfl998h . 

For Langmuir oscillations, Landau damping in a hot, isotropic, unmagnetized plasma 
is controlled by the spectral index while for ion-acousti c waves Landau dam ping is more 
sensitive to the ion temperature than the spectral index dOureshi et al. 1120061) . Thus, Lan- 
dau damping growth rates of l ong wavelengths Langmuir m odes become much larger for 
a Lorentzian (Kappa) plasma dThorne and Summers! 1 1 99 lh limiting the existence of the 
(weakly damped) Langmuir waves in space plasmas with Kappa distributions, to a narrow- 
band just above the electron plasma frequency. The significant increase of spatial Landau 
damping by small K electrons is also demonstrated for spatially damped pla sma waves gen - 
erated by a planar grid electrode with an applied time harmonic potential dPodestal 12005). 
Hybrid models of Maxwellian plasmas partially populated by hot fc-components can pro- 
vide better fits to the observations and experiments than the simple uniform Maxwellian 
dHellberg et al.ll2000d : lMace et alJll999l) . 

The relativistic effects on dispersion and Landau dam ping of Langmu ir waves in a rel- 
ativistic Kappa-distributed plasma have been studied by IPodestal d2008l) . The relativistic 
dispersion relations derived have been used to compute the damping rates and phase speeds 
for plasma waves in the solar wind near the Earth orbit. It was found a good match for the 
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electron velocities in the superhalo with the phase speed of weakly damped plasma waves, 
and thus providing a plausible mechanism for their acceleration. 

The e xistence conditions and ch a racteristics of ion- aco ustic solitary wave s have been 
studied bv lAbbasi and Paiouhlj2008h ; ISaini etaD j2009h and lChuang and Hau I d2009h show- 
ing that Kappa-distributed electrons are not favorable to the existence of these solitons. A 
comparative study of Langmuir waves, dust ion acoustic wa ves, and dust - acous tic waves 
in Maxwellian and Kappa distributed plasmas is presented in IZaheer et al.l d2004h . Landau 
damping rate of dust acoustic waves in a dusty plasma modeled by a Kappa distribution for 
electrons and ions and by a Maxwellian for the dust grains has bee n found to be dependent 
on the spectral index K as well as the ratio of ion density to electron (lLeel2007h . Dust acous- 
tic solitons have also been studied in pla smas with Kappa-distributed electrons and/or ions 
and c old negative or positive dust grains (Bal uku and Hellberg 2008; Younsi and Trib eche I 
l2008h or with non thermal ions having Kappa-vortex-like velocity distributions functions 
feamel et al1l2009h . 

5.3.2 Magnetized plasmas 

The general dielectric tensor for magnetoplas mas comprising compon ents with generalized 
Lorentzian distributions has been calculated bv lSummers et al.l i 19941) for arbitrary oriented 
wave- vectors. Applying to the electrostatic or electromagnetic waves propagating parallel to 
the ambient magnetic field, simple dispersion r elations can be deri ved dSummers and Thornd 
1 19911 ; IXue et al. Ifl993l ; ISummers et all 1 1994 lMacdll996l. 1 19981) in terms of the modified 
plasma dispersion function Q). 

This dielectric tensor has also been simplified for a 3-dimension al isotropic Kappa dis- 
tribution in a form similar to that obtained by Trubnikov (Mace 1996), and for a Kappa loss- 
cone distribution with applications to a large variety of space plasmas like th e solar wind , 
magnetosheath, ring current plasma, and the m agnetosp heres of other planets (Ma cell 19961; 
Xiao et al. ll998l ; |Pokhotelov et al. l2002l ; lxiao l2006l ; lxiao et al. l2006al ; lsinghal and Tripathi I 
2007h . 

For wave-vectors oblique to the magnetic field, iMacd d2003l , 12004]) has described the 
generalized electron Bernstein modes in a plasma with an isotropic Kappa velocity distribu- 
tion. In a hybrid Kappa-Maxwellian plasma, unlike the uniform Maxwellian plasma, the dis- 
persion properties of the oblique electromagnetic waves were found to be markedly changed 
from an elaborate study including effects of the Kap pa value, the propag ation angle, and the 
temperature anisotropy on dispersion and damping dCattaert et al. 1 12007). Notice that, be- 
cause of the anisotropy of the countours in the velocity space, such a Kappa-Maxwellian 
distribution is unstable in an overdense plasma near the electron-cyclotron frequency even 
when the parallel and perpendicular temperatures are equal. 

The dispersion relations for low- frequency hyd romagnetic waves in a Kappa distributed 
plasma has recently been derived by I Basil I (|2009) showing that both the Landau damping 
and the transit-time damping (magnetic analogue of Landau damping) of the waves are 
enhanced in the suprathermal region of the velocity space. 

5.4 Anisotropic Kappa distributions 

The effects of anisotropic Kappa distributions have also been investigated since the anisotropic 
velocity distributions are most probably at the origin of nonthermal emissions in astrophys- 
ical sources, and the magnetic field fluctuations in space plasma. 
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Fig. 4 Dependence of dispersion curves (panel a) and growth rates (b) for the electron cyclotron instability 
(R mode) on the spectral index K — 2 (red dotted lines), 4 (blue dashed lines), and K — > °° for Maxwellian 
plasmas (black solid lines). These are exact numerical solutions for parameters intended to simulate (low 
altitude) solar wind conditions: plasma temperature T e = 2.5 x 10 7 K, a temperature anisotropy T±/T e ~ 4, 
and gyrofrequenc y \Q e \ = 0.1a) pe . T he coordinates are scaled as Wr= w r / co pe Wi= to,-/ to pf and K = kc/w pe . 
(Adapted from lLazar et aI1fe008alU 



When both electrons and ions are modeled by an anisotropic distribution of bi-Kappa 
type, I Summers and Thornd Jl99lh have derived the general dispersion relation for the par- 
allel electromagnetic modes, right-handed (R mode) and left-handed (L mode) circularly 
polarized, in terms of the modified plasma dispersion function The effect of suprather- 
mal particles on the stability of these modes largely varies depending on the shape of the 
distribution function, and the mode frequency, whether it fits to thermal Doppler shift of the 
electron gyrofrequency (the whistler instability driven by the cyclotron resonance with elec- 
trons) or the ion gyrofrequency (the electromagnetic ion cyclotron instability). Thus, while 
the growth rates of the electron cyclotron instability (R-mode lower branch, CO, < \£2 e \) 
become lo wer in a bi-Kappa p lasma than for a bi-Maxwellian with the same temperature 
anisotropy jLazar et al, 2008a), also see Fig. [4] b, at smaller frequencies (c o r <C \£i e \) . the 
whistl er growth rates become higher dMacel 1 998l : lTripathi and Singhall2 008). Cat taert et al. I 
(120071) have also shown that, unlike a bi-Maxwellian plasma, the low-frequency whistler 
modes in a Maxwellian-Kappa plasma (described above), can be stable to the temperature 
anisotropy. Their study includes effects of varying fc for both underdense and overdense 
plasmas, and for both parallel and oblique propagation. 

The loss cone bi-Lorentzian distribution, which allows to plasma populations to have 
anisotropic temperatures and a loss cone, has been used extensively (although the latter 
is largely inconsequential in mod els for wave propagation parallel to the magnetic field) in 
space and laboratory applicat i ons i Macdl998:lTripathi and Singhal 2008l lTripathi and Misral 



2000(). ISinghal and Tripathi ] J2007n obtained the components of the dielectric tensor for 
such a distribution function and made parametric studies of the effect of fc-index, the loss- 
cone index (which is , in g eneral different from fc), and different temperature anisotropies 
( iTripathi and Singhall l2008). Whistler mode instability in a Lorentzian (Kappa) magneto- 
plasm a in the presence of perpe ndicular AC electric field and cold plasma injection was stud- 
ied bv lTripathi and Misral (2000). An unperturbed Lorentzian distribution has also been used 
for studying the effect of a cold plasma beam on the electromagnetic whistler wave in the 
prese nce of a perpendicular AC electric field in the Earth's atmosphere dPandev and Pandevl 
120081) . 
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The Kappa-loss-cone (KLC) distribution function obeys a power-law not only at the 
lower energies but al s o at th e relativistic energies. A relativistic KLC distribution has been 
introduced by Ixiao I ( 120061) for an appropriate characterization of the energetic particles 
found in planetary magnetospheres and other plasmas, where mirror geometries occur, i.e., 
a pronounced high energy tail and an anisotro py. The field - aligned whis tler growing modes 
in space plasmas have been investigated by Ixiao et al. I |2006a) and Izhou et all ([2009) 
applying relativistic treatments for relativistic Kappa or KLC distributions. The threshold 
conditions f or the w histler instability in a Kappa distributed plasma have been derived by 
Xiao et al. 1 d2006bh . Numerical calculations were carried out for a direct comparison be- 
tween a KLC distribution and the current Kappa distrib ution. The KLC wa s also adopted to 
model the observed spectra of solar energetic protons dXiao et al~ll2008bh . Recent studies 
dXiao et al. ll"2008cJ : Izhou et al. 112009b have introduced a generalized relativistic Kappa dis- 
tribution which incorporates either temperature anisotropy or both loss cone and temperature 
anisotropy. 

Purely growing (co ~ 0) mirror modes and low frequency electromagnetic ion cyclotron 
waves are widely detected in the solar wind and magnetosheath plasmas being driven by an 
ion temperature anisotropy, T± > 7ji (or pressure, px > P\\)- The effects of the suprathermal 
tails on the threshold conditions and the linear growth rates of these instabilities have widely 
been investigated in the last d ecade dLeubner and Schup fer 2000, 200ll . l2002hlGedarin et al.l 
l200ll : |Pokhotelov et al. l2002n . An universal mirror wave-mode threshold condition for non- 
thermal space plasma environments was obtained by iLeubner an d Schupfer d2000L 1200 lL 
2002). The li near theory of th e mirro r instability in non-Maxwellian space plasmas was de- 
veloped by IPokhotelov et al. I d2002h for a large class of Kappa to suprathermal loss cone 
distributions in view of application to a variety of space plasmas like the solar wind, mag- 
netosheath, ring current plasma, and the magnetospheres of other planets. Thus, while tran- 
sition to nonthermal features provides a strong source for the generation of mirror wa ve 
mode activity, reducing drastically the instability threshold ( lLeubner"an d Schupfer 20q2), a 
more realistic presence of suprathermal tails exclusively along the magnetic field (parallel fc- 
distribution) counteracts the growth of the mirror instability and contributes to stabilization 
jPokhotelov et al. I [2002 ) . In the nonline ar regime, solitary struc tures of the mirror waves 
occur with the shape of magnetic holes jPokhotelov et al. H OOPS') suggesting that the main 
nonlinear mechanism responsible for mirror instability saturation might be the magnetic 
trapping of plasma particles. 

The ion cyclotron wave instability driven by a temperature anisotropy (7j_/Iji > 1) of 
suprathermal ion s (protons) modeled with a typical Kappa distribution is investigated by 
IXue et al. I ill 993b : Ixiao et all d2007ah : IXueet al. I dl996al lbh for solar wind co nditions and 
for mag netosphere. The threshold condition for this instability is determined bv lXiao et al. I 
(2007b). As the spectral index K for protons increases, the maximum growth rates of R, L 
modes decrease ilSummers and Thornelll992l : IXue et al. ll993l : lDasso et al.l2003l : IXue et al71 
1996b) and the instability threshold generally decreases and tends to the lowest limiting 
values of the bi-Maxwellian (fc — > °°) dXiao et al. Il2007bh . The corresponding enhancement 
in the growth rate of L mode waves in planetary magnetospheres is less dramatic, but the 
Kappa distribution tends to produce a significant wave amplification over a bro ader range 
of fre quency than a Maxwellian distribution with comparable bulk properties ( IXue et al~l 
1 19931) . The damping and growth rates of oblique waves are also lower for K appa distri- 
bution s, but differences become less important for nearly perpendicular waves llXue et al. I 
Il996ah . 

In the opposite case, a surplus of parallel kinetic energy, T±_ < 7ji (or pressure p± < pu) 
will excite other two kinetic instabilities: the firehose instability propagating parallel to the 
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Fig. 5 Dependence of dispersion curves (a) and growth rates (b) for firehose instability on the spectral index 
K = 4 (red dotted lines), 5 (blue dashed lines), and K — > °° for Maxwellian plasmas (black solid lines). These 
are numerical exact solutions obtained with parameters estimated for solar flares: plasma density n = 5 x 10 10 
cm~ 3 , electron and proton temperature T e ± = T p ± = T p u ~ 10 7 K, the electron temperature anisotropy 
T e t ±/T e II ~ 20 (in a), and _L and || denote directions with respect to the ambient magnetic field Bp = 100 G. 
The co ordinates are scaled as Wr= (Q r /Q p , Wi= Q)j/Q p and K = kc/w p . (Adapted from Lazar and Poedts 
<2009h ) 
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Fig. 6 Dependence of growth rates for Weibel instability (panel a), filamentation (b) and two stream insta- 
bility (c), on the spectral index K — 2 (red dotted lines), 4 (blue dashed lines), and K— > °° for Maxwellian 
plasmas (black solid lines). Here the parameters are typical for solar wind conditions: plasma temperature 
r c ~2x 10 6 K, a temperature anisotropy T±/T e ~ 4 (in a), and symmetric counterstreams in b and c, with the 
same temperature and the bulk velocity vo = 0.1c (where c = 3 x 10 s m/s is the speed of light in vacuum). 
The c oordinates are scal ed as Wi=ffl,/co / , e and K = kc/w pe in panels a and b, and K = k\'o / co pe in c. (Adapted 
from lLazar et al. I <2008bl3) .) 
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magnetic field lines and with maximum growth rates of the order of ion gyrofrequency 
dLazar and Poedtsll2009h . and the Weibel like instability propagating perpendicular to the 
hotter direction (in this case also perpendicular to the magnetic field lines), and which, 
in general, is much f aster reaching max i mum g rowth rates of the ord er of electron (or 
ion) p lasma frequency dLazar et al. I2008bll200 9a). Recent investigations dLazar and Poedts 1 
2009) of the electron firehose instability driven by an anisotropic electron distribution of bi- 
Kappa type, have proven that, by comparison to a bi-Mawellian, the threshold increases, 
the maximum growth rates are slightly diminished and the instability extends to large wave- 
numbers (see Fig. |5J. Instead, a more important reduction has been found for the growth 
rates of the electron Weibel instability (Fig. |6] a) in a non-magnetized or weakly magn e- 
tized plasma with bi-Kappa distributions dZaheer and Murtaz a 2007; Lazar et al . 2008b, a, 
l2009ci) . 

Despite the similar features of the electromagnetic filamentation instability, which is a 
Weibel-like instability that grows in a counterstreaming plasma or a beam-plasma system 
perpendicular to the streaming direction, the effect o f suprathermal populations is opposite 
enhancing the filamentation growth rates (Fig. [6] b) dLazar et al. Il2008bl . l200 9ar). Extended 
investigations have included counterstreaming plasmas with an internal bi-Kappa distribu- 
tion, when the filamentation and Weibel effects cumulate leading again to increased growth 
rates but only for plasmas hotter in the streaming direction. Otherwise, if counterstreaming 
plasmas are hotter in perpendicular direction, t he effective anisotropy dec reases, diminish- 
ing the growth rates of filamentation instability dLazar et al. I2008bl .l2009a). However, in this 
case, two other unstable modes are expected to arise, both along the streams: a Weibel-like 
electromagnetic instability and a two-stream electrostatic instability, which is, in general, 
faster than Weibel. Furthermore, suprathermal populations enhance the electrostatic insta- 
bility leading to larger growth rates for lower K (Fig. [6] c). The same behavior has been 
observed for the modified two-stream instability driven by the relative motion of ions, as- 
sumed Kappa-distributed, with respect to the el ectrons (assumed Maxw ellian) : maximum 
growth rates decrease with K and with plasma /3 (Langma vr et al. 1 12005). 

Recently, iBasu (2008, 2009) provided a systematic study for the stability of a magne- 
tized plasma at low frequencies and in various limits of low or high plasma beta (plasma 
pressure/magnetic pressure), showing that the threshold values for the excitation of the un- 
stable hydromagnetic waves in Kappa distribution plasma are increased as a consequence of 
enhancing the resonant wave-particle damping. 



6 Outlines and perspectives of Kappa distributions 

The Kappa function has been proved to be a convenient tool to describe plasma systems out 
of thermodynamic equilibrium. Since the fast particles are nearly collisionless in space plas- 
mas, they are easily accelerated and tend to produce nonequilibrium velocity distributions 
functions with suprathermal tails decreasing as a power law of the velocity. Thus, models 
based on the kappa distribution allow to analyze the effects of the suprathermal particles 
and to fit distributions for ions and electrons measured in situ, in the solar wind and in the 
magnetosphere of the planets. 

Major consequences follow the presence of these suprathermal particles, and especially 
the velocity filtration which makes the kinetic temperature increase upwards. This mech- 
anism has been proposed to explain the heating of the corona. The suprathermal particles 
increase also the escape flux in planetary and stellar wind and can explain the accelera- 
tion of the fast solar wind. For low values of fc, the heat flux changes of sign compared to 
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Spitzer-Harm so that heat can flow from cold to hot. The future solar missions (Solar Orbiter 
and Solar Probe) should improve the observations concerning the presence of suprathermal 
particles in the corona. 

Valuable theories have been proposed concerning the origin and the fundamental phys- 
ical arguments for a kappa family distribution function. The universal character of these 
distribution functions suggests that they can be attributed to a particular thermodynamic 
equilibrium state related to the long range properties of Coulomb collisions. Thus, Kappa 
distributions generalizes the notion of equilibrium for collisionless plasmas far from thermal 
(Boltzmann-Maxwell) equilibrium, but containing fully developed quasistationary turbulent 
fields. A new statistical mechanical theory has been proposed extending Gibbsian theory 
and relaxing the independence of subsystems (no binary collisions) through introducing a 
generalized Kappa function dependence on entropy. In the absence of binary collisions, the 
Kappa parameter has been found to be a measure for the strength of the subsystem cor- 
relations introduced by the turbulent field fluctuations. Particular forms for a nonextensive 
(superadditive) entropy and for temperature have been derived to satisfy the fundamental 
thermodynamic relations and recover exactly the family of Kappa distributions from 0. 

An isotropic Kappa distribution will therefore be stable against the excitation of plasma 
instabilities. This function replaces the Boltzmann-Maxwell distribution in correlated col- 
lisionless equilibria of plasma particles and turbulent fields. But such dilute plasmas eas- 
ily develop flows and temperature anisotropies making Kappa distributions to deviate from 
isotropy and drive kinetic instabilities. In space plasmas embeded by the interplanetary mag- 
netic field, heating and instability are, in general, resonant, and both will be enhanced in 
Kappa distributed plasmas because there are more particles available at high energies to res- 
onate with waves. Revisiting transport theories and calculation of the transport coefficients 
in solar environments is therefore an important task, and a significant progress is expected 
to be done involving Kappa distributions. 

Many solar models have been developed on a questionable existence of a Maxwellian 
equilibrium, so that, now, the presence of Kappa distributions asks for more realistic rein- 
terpretations that are expected to provide better fits to the observations. The new techniques 
developed for measuring the electron density, temperature, and the suprathermal parameters 
will also offer important clues to understanding transport properties in space plasmas. 
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